Self-energy Effect of Superconducting Energy Gap in Point-Contact Spectra of MgB 2 
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In strong-coupling superconductors with a short electron mean free path the self-energy effects in 
the superconducting order parameter play a major role in the phonon manifestation of the point- 
contact spectra at the above-gap energies. We derive asymptotic expressions of the phonon structure 
in MgB2 in the case eV 3> A for tunnel, ballistic, and diffusive point-contacts and show that these 
expressions not only qualitatively, but also semi-quantitatively correspond to the measurements of 
the phonon structure in the point-contact spectra for the 7r-band of MgB2 c-axis oriented thin films. 
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Introduction. 

It is commonly accepted that the mechanism of super- 
conductivity in recently discovered MgB? [I| is due to 
the electron-phonon interaction (EPI) [2J. Among the 
known s — p metals, the record breaking critical temper- 
ature (T c « 40 K) and the unusual two-band character in 
its electronic structure make this compound very inter- 
esting for detailed study. MgB2 crystallizes in a hexag- 
onal lattice with alternating planes of Mg and B. Up to 
now, there has been no detailed experimental determi- 
nation of its EPI spectral function. The reason is that 
EPI in MgB2 is strongly anisotropic, and single crystals, 
well suited for anisotropic measurements by means of 
tunneling or point-contact spectroscopy, are still under 
development. As for the tunneling spectroscopy, in or- 
der to extract the so-called "driving force" for high T c in 
this compound caused by EPI along the a6-plane [3j , one 
needs a plane tunnel junction oriented perpendicular to 
this direction, which is not an easy task. On the other 
hand, the tunneling spectra in the c— direction are sup- 
posed to be so weak (see below) , that the precise standard 
procedure for extracting the EPI function is difficult 
to apply. In Ref. [j| it was shown that the inversion of 
the Eliashberg equation for a multiband superconductor 
is a mathematically ill-defined problem. The previous 
attempt to obtain a solution on polycrystalline samples 
contains an uncontrolled mixture of the two bands with 
a predominant contribution of the 7r-band Q . The inter- 
action of the two bands to a great extent determines the 
corresponding EPI functions. Their intensity depends on 
electron mean free path in each band separately. On the 
one hand, in some experiments the two bands can be 
considered as being independent of each other, but, on 
the other hand, both of them have the same T c , despite 
quite different EPI. The layered structure and strongly 



anisotropic EPI allow us to hope that the peculiarities in 
MgB2 may shed some light on the mechanism of high-T c 
superconductivity in copper oxides. 

In the previous paper pj it was shown qualitatively 
that the phonon singularities in the point-contact spec- 
tra of the c-axis oriented MgB2 thin films were due to 
the non-linear dependence of excess current I exc (eV) at 
energies higher than the superconducting energy gap A. 
In the present contribution we derive asymptotic expres- 
sions for differential conductance of N — c — S (c stands 
for " constriction" ) point-contacts and show that quanti- 
tatively (within a factor of order of unity) the calculated 
point-contact spectra agree with the experimental data 
for the c-axis orientation of MgB2 , where only the 7r-band 
is visible. 

Theory. 

The electronic structure of MgB2 consists of two 
groups of bands 0, Q : a P a i r °f approximately isotropic 
3D 7r-bands, and a pair of strongly anisotropic 2D er- 
bands whose characteristics are measurable for the point 
contact oriented within a few degrees around the a6-plane 
H . As it was argued in Q , the variation of the supercon- 
ducting gap inside the a- or the 7r-bands can hardly be 
observed in real samples due to intraband impurity scat- 
tering. Therefore, superconducting properties of MgB 2 
can be described by an effective two-band model, where 
each group of the a- and the 7r-bands is characterized 
by the corresponding order parameters A^. and A T . The 
gap functions A CTjTr in MgB2 in the Matsubara representa- 
tion, suitable for the description of thermodynamic prop- 
erties, were calculated in Ref. 01 from the solution of 
the Eliashberg equations using the theoretically calcu- 
lated EPI spectral functions for the effective two-band 
model. Here we extend this approach to obtain the com- 
plex gap parameters A (7)ff (£') in MgB2 as a functions of 
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FIG. 1: The self-energy effects in the complex order param- 
eter of superconducting MgF^. The upper and lower panels 
are for the it- and the cr-band, respectively. 



the real energy E. The results are shown in Fig. ^ an d 
demonstrate the self-energy effects important for the de- 
scription of the transport properties across MgB2 point 
contacts. 

If the point contact is clean, then the electron mean 
free path / is greater than the size of the constriction d, 
and the current flows in the ballistic regime. The differ- 
ential conductance of S — c — N contact at T = obeys 
the following expression 
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eV, 



where Rn is the normal state resistance. 

In the dirty contact (I <C d), the regime of current 
flow is diffusive and the corresponding expression for the 
differential conductance follows from formula (21) from 
Ref. [13 for the current- voltage dependence: 
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It is well known that in the tunnel regime the differ- 
ential conductance is proportional to the quasiparticle 
density of states jj| : 



A 2 (e)' 



eV. (4) 



We assume T = in these formulae, but they can be 
generalized to a finite temperature as well. 

These formulae acquire a very simple form in the limit 
eV 3> A. For tunnel, ballistic and diffuse point contacts 
they are as follows: 
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eV, respectively. Thus, the dl/dV- 
characteristics of any point contact (tunnel, ballistic, dif- 
fusive) are roughly similar below 80 meV, where the con- 
tribution of ImA(e) is negligible. 

In order to compare the theoretical predictions with 
the experimental spectra we express the calculated spec- 
tra as the second harmonic (i.e., rms voltage, V2(V)) 
of the small first harmonic modulation voltage V\ by 
expanding the V(I) characteristic in the Taylor series. 
Namely, we express the V2-signal as 



V? d ( dl\ 

^dV [ RN dv) 



(8) 



assuming that changes in differential resistance is small 
and, correspondingly, the modulation voltage is fixed at 
V\. Any other spectra can be reduced to the particular 
value of V\ using the equation JSJl. In what follows we 
choose V\ = 3 meV, which is close to the value used in 
the experiments. 

In Fig. |21(a), (b) the calculated by formulae JSJl, I©, 
0, and (JSJl second harmonic voltages V^(Vi = 3 meV) 
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FIG. 2: The calculated second harmonic signal V2(V) propor- 
tional to the second derivative of the / — ^-characteristic of 
the N—c—S point-contact with a ballistic or a diffusive regime 
of current flow (see text). For comparison, the corresponding 
curve for the tunnel junction N — I — S is shown together 
with the BCS tunneling density of states. Calibration of V2 
is reduced to modulation voltage Vi=3 mV. 



are shown for the 7r— and a— bands of MgB2. These ex- 
treme regimes of current flow: tunnel, ballistic, and dif- 
fusive, are calculated using A(eV) functions from Fig. Q 
One can see a small amplitude (of the order of 1 uV) 
of the phonon structure for the n— band, and a much 
stronger E 2g phonon mode singularity at about 80 meV 
for the cr-band. Up to eV ~ 80 meV the self-energy 
structure is determined by ReA(e). The amplitude of 
phonon structure is different, reflecting the correspond- 
ing prefactors 1/2 — > 1/3 — > 1/4 in the row of current 
flow regimes: tunnel — *■ diffusive — > ballistic. On the 
other hand, large differences in shape are observed in the 
energy range 80-^120 meV, where ImA(e) becomes appre- 
ciable. We shall see below that just at the high energy 
edge the experimental curves show some variation due to 
the uncontrollable changes of scattering in the contact 
region. 

Comparison with experiment. 

Since the theoretical spectra are odd with regard to 
the changes of the dc voltage sign, we extract the odd 
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FIG. 3: The extraction of an odd part of the experimental 
point-contact spectra shown for the middle curve in Fig. 4 
(a). The thinner curve shows the raw data shifted vertically 
down for clarity. 



part of the raw data as shown in Fig. [3J 



V 2 ° M = ~[V 2 (eV)-V 2 (-eV)] 



(9) 



In Fig. 2] we compare 3 odd parts of the experimental 
spectra for different point-contacts which are perpendic- 
ular to the surface of the c-oriented thin film [13j . The su- 
perconducting energy gaps An measured by the Andreev 
reflection spectroscopy from <iF/<iZ-characteristic lie in 
the range 2.1-^2.6 meV. These (iV/dl-curves show almost 
no traces connected with the large gap at about 7 meV. 
Hence, for this orientation we probe only the 7r-band ||. 
We compare the experimental curves (Fig. 0] (a)) with 
the 7r-band calculated spectra (Fig. 0| (b)), since in the 
experiments the regime of current flow is among the ex- 
tremes calculated by the theory. One can see that not 
only the shape of the experimental spectra corresponds 
well to the theoretical ones (with minor deviations), but 
also the amplitude of the structure has a proper order of 
magnitude (~ 0.1 % of R N ). From the dV/dI(V)Andreev 
reflection spectra (see the caption in Fig. ??) one can 
suggest that our contacts are closer to the tunnel regime 
in the series of the experimental curves with resistances 
49 — * 36 — > 80 Q. To make the estimation more quanti- 
tative, one needs an interpolation formula similar to the 
finite barrier parameter Z in the BTK-theory [hif. 

The question arises as to why we did not observe a 
much stronger self-energy S-type-structure in Ref. |l5l |. 
where the point-contact spectrum of MgB2 single crystal 
was measured along et6-planc? The answer is that in that 
case the nonequilibrium phonon generation is so strong 
that superconductivity in the contact region is destroyed, 
and correspondingly, the excess current dramatically de- 
creases around the E 2g phonon mode energy. This leads 
to a large maximum in the dV/dI(V) characteristic and 
the N-shape singularity in V 2 (V). The maximum can be 
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FIG. 4: Comparison of three different experimental point- 
contact spectra (panel (a)) in the c-direction (Ref. 0) with 
the theoretical ones (panel (&)) calculated by asymptotic for- 
mulae f |5|. ©, 0, and based on superconducting or- 
der parameter A(e) from Fig. Q All the curves are dis- 
placed vertically for clarity and reduced to modulation voltage 
Vi = 3 mV. In panel (a), the normal state resistances 49 and 
36 O correspond to zero-bias Ro= 45 and 43 Q, respectively, 
being the same as in Fig. 2 (b) of Ref. 0. For R N = 80 Q 
7?o=150 Q. The experimental curves are smoothed by several 
mV not disturbing the phonon structure. T—4.2 K, H=0. 

easily seen in the Ui(U)-characteristic (note, for example, 
the inset in Fig. 1 of Ref. [l^)- Normally, this destruc- 
tion of superconductivity masks the self-energy structure 
effectively It can be identified by the strong dependence 
of energy position on external parameters: magnetic field 
and temperature. 
Summary. 

In conclusion, we have derived the simple asymptotic 
formulae for the order parameter self-energy effects in 
the superconductin g po int-contact. They can be used in 
standard programs Il7| to solve the Eliashberg equa- 
tions [4j for quantitatively derivation of electron-phonon- 
interaction spectral function providing the structure of 
point-contact is established. In MgB2, we applied them 



in c-direction for obtaining point-contact spectra in ir- 
band. The close similarity between the calculated and 
measured point-contact spectra in the c-direction mani- 
fests the validity of the calculated EPI spectral function 
in the 7r-band in Ref. [Icj . 
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